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A SIMPLIFIED METHOD FOR CALCULATING LAMINAR HEAT TRANSFER 
OVER BODIES AT AN ANGLE OF ATTACK 

By Fred R. DeJarnette and Ruby M. Davis 
Langley Research Center 

SUMMARY 

A simplified method is developed for calculating the ratio of the local to stagnation- 
point heat-transfer rate for bodies at an angle of attack with only the free-stream Mach 
number and the ratio of specific heats required as inputs. The viscous problem is sim- 
plified by using the axisymmetric analog for three-dimensional boundary layers (small 
cross flow in boundary layer) in conjunction with Lees’ laminar heating rate for axisym- 
metric bodies. An approximate technique for determining the geometry of the inviscid 
surface streamline is presented. In this technique the direction of a streamline at a 
given point on the body is taken as the direction of the resultant of the free-stream veloc- 
ity vector minus its normal component. For determining the inviscid surface proper- 
ties, the modified Newtonian pressure distribution is used with isentropic flow along the 
surface. 

The heat transfer over spherically blunted cones with 15° and 30° semiapex angles 
at angles of attack of 0°, 10°, and 20° and a free-stream Mach number of 10.6 was com- 
puted by the present method. The values obtained agreed well with experimental data. 

INTRODUCTION 

Theoretical methods for predicting the aerodynamic heat transfer to symmetric 
bodies at an angle of attack and to lifting bodies are required for the proper design of 
reentry heat-protection systems. A review of the literature (see, for instance, ref. 1) 
shows that previous methods are restricted to small angles of attack or to yawed infinite 
cylinders, or that they involve complicated numerical techniques to solve the governing 
flow equations. This paper presents a relatively simple method for computing the laminar 
convective heating rates over basic hypersonic configurations at an angle of attack and 
with a highly cooled surface. Radiative heat transfer is not considered here, nor is the 
effect of ablation gases on the convective heating rates. 

The complexity of the three-dimensional differential equations governing the invis- 
cid as well as the viscous flow field over a body at an angle of attack makes the use of 
simplifying approximations desirable so that tractable solutions may be obtained. A 
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substantial simplification to the viscous flow field equations is achieved through the axi- 
symmetric analog for three-dimensional boundary layers (refs. 2 and 3). This axisym- 
metric analog transforms the general three-dimensional boundary-layer equations into 
the same form as those for an axisymmetric body at an angle of attack of 0° when the 
coordinate directions are taken along and perpendicular to the inviscid body streamlines, 
and when the cross flow (i.e., the component of boundary- layer flow normal to an inviscid 
streamline and along the body surface) is assumed to be small. Therefore, if the pres- 
sure distribution and streamline geometry are known, the heat transfer may be calculated 
along a streamline by any method applicable to a body of revolution at an angle of attack 
of 0°; such as the method of Lester Lees (ref. 4). 

The major difficulty in applying the axisymmetric analog to angle -of- attack problems 
is that the inviscid solution on the body surface (surface pressures, velocity, and stream- 
line geometry) is required. Ordinarily, the inviscid solution would have to be obtained 
by some numerical technique, which is a major undertaking in itself for bodies at an 
angle of attack. Therefore, some approximations in the inviscid solution are, in many 
cases, necessary. 

In the present method, the direction of an inviscid streamline at a given point on the 
body is obtained by assuming that the free-stream velocity vector loses its normal com- 
ponent upon striking the body. Thus, the direction of a streamline is taken as the direc- 
tion of the resultant of the free-stream velocity vector minus its normal component at a 
given body point. Once the expression for the streamline direction is obtained, a first- 
order ordinary differential equation for the geometry of the streamline may be written. 
Then, another ordinary differential equation for the scale factor (corresponding to the 
body radius in the axisymmetric analog and often referred to as the streamline diver- 
gence function) is obtained. Both these differential equations may be readily integrated 
in closed form for some simple body shapes; and, for more complicated body shapes, 
they may be integrated by a simple numerical scheme such as the Runge-Kutta method. 

The magnitude of the flow velocity along a streamline is determined directly from 
the pressure distribution under the assumption that the inviscid flow on the surface is 
isentropic. In the present paper, modified Newtonian pressure distributions are assumed 
for convenience. 

The accuracy of the present method was partially assessed by computing the ratio 
of local to stagnation-point heat-transfer rate on spherically blunted 15° and 30° cones 
at a free-stream Mach number of 10.6, a ratio of free-stream specific heats of 1.4, and 
angles of attack of 0°, 10°, and 20°. The results are shown to compare well with experi- 
mental data. 
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SYMBOLS 


A function defined by equation (10) 

e n ,es unit vectors normal (outer) to the body surface and along inviscid surface 

streamline, respectively 

e^ unit vector normal to e n and e s , defined by equation (12) 

e r ,e^ unit vectors in r and <f) directions, respectively 

G function defined by equation (3) 

H scale factor for £ coordinate 

h scale factor for /3 coordinate 

i,j,k unit vectors in x, y, and z directions, respectively 

k slope of cone surface 

l metric defined by equation (5) 

Moo free- stream Mach number 

n coordinate normal to body surface 

p local static pressure 

p Q stagnation pressure behind normal shock wave 

q w local heat-transfer rate 

R nose radius 

r body radius 

rb body base radius 



s distance along a streamline, measured from the stagnation point 

u e inviscid velocity along body-surface streamline 

Voo undisturbed free-stream velocity 

Vco undisturbed free-stream velocity vector 

x,y,z Cartesian coordinates (fig. 1) 

a angle of attack 

/3 coordinate measured along body surface, perpendicular to inviscid streamlines 

y^ ratio of specific heats in free stream 

y mean ratio of specific heats behind bow shock wave 

6 inclination of body surface with respect to body axis 

£ coordinate measured along a streamline, H d£ = ds 

4> angular coordinate (fig. 1) 

angle between Voo and -e n , defined by equation (20) 

Subscripts: 

i sphere- cone interface 

o stagnation point 

x indicates that derivative is taken with x held constant 

/3 indicates that derivative is taken with (3 held constant 

A prime indicates a derivative with respect to x. 
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ANALYSIS 


Problem Description 


For axisymmetric blunt bodies at an angle of attack of 0° and with a highly cooled 
surface, the ratio of the local heat-transfer rate q w to the stagnation-point heat-transfer 
rate q w Q given in the notation of the present report, was derived in reference 4: 


*lw,o 


(P/P 0 )(V V oo ) rRl/2 

(P/P 0 )(u e /Voo)r 2 ds 


2G 


( 1 ) 


The modified Newtonian pressure distribution combined with the assumption of isentropic 
flow along the body surface yields the following results (ref. 4): 


P 


P 


o 




cos +. 



( 2 ) 


G = 




2 \ 

(roo - l ^ 2 / 1 




( 3 ) 


u e J 

1 1 2 

/ p \( y ~ y\ 


(roo ~ l)M,o 2 

' p ° Jj 


(4) 


According to the axisymmetric analog (refs. 2 and 3), equations (1) to (4) are also 
applicable to any inviscid surface streamline on an asymmetric body or an axisymmetric 
body at an angle of attack if s is the distance measured along the streamline and r is 
replaced by the scale factor h corresponding to the coordinate 0 measured along the 
body surface and perpendicular to the streamline. The assumption of zero cross flow, 
which is implied in the axisymmetric analog, was shown in reference 5 to be applicable 
to highly cooled blunt bodies in hypersonic flow. 


The present analysis will consider a spherically blunted axisymmetric body at an 
angle of attack. The basic method presented here is applicable, however, to any three- 
dimensional body at an angle of attack. The coordinates £, j3, and n (when n is 

measured from the surface along a straight line normal to the surface) form an ortho- 
gonal coordinate system whose metric is 


dZ 2 = H 2 d£ 2 + h 2 dj3 2 + dn 2 


(5) 


where H d| = ds. 
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On the body surface n = 0, H = H(x,0), and h = h(x,0), where x is the distance along 
the body axis and 0 is the circumferential position. (See fig. 1.) The coordinate 0 
is constant along a given inviscid surface streamline. 

If the modified Newtonian pressure distribution (eq. (2)) is accepted as being suffi- 
ciently accurate, then the heating- rate ratio for the body at an angle of attack may be 
computed from equations (1) to (4) along any inviscid surface streamline, provided the 
streamline geometry and scale factor h are known (H is not required). 


Streamline Geometry 

In order to determine the streamline geometry and scale factor h, it is assumed 
here that the direction of an inviscid surface streamline may be taken as the direction of 
the resultant of the free-stream velocity vector minus its normal component at a given 
point on the body surface. With e s defined as a unit vector in the direction of a stream- 
line, this assumption gives 

g _ Voo ~ ( Voo • e n) e n 
|v oo -(v oo .e n )e n | 

where e n is a unit vector normal (outer) to the surface, which is given by 

•ir 1 + j cos 0 + k sin 0 


(6) 


e n - 


(1 + r'2) 1 / 2 


(7) 


The quantities i, j, and k are unit vectors in the x, y, and z directions, respec- 
tively (fig. 1). As shown in figure 1, lies in the xy- plane; hence 

Voo = V DO (i cos a + j sin a) 


( 8 ) 


where a is the angle of attack. Equations (7) and (8) may be utilized in equation (6) to 
obtain 


e s = 


(i + r'e r )(cos a + r' sin a cos 0) - e^( 1 + r'2)sin a sin 0 


(1 + r'2)V2 A 


where 


A = jsin20 + (r' sin a + cos a cos 0)^j 


1/2 


(9) 


( 10 ) 


and e r and e^ are unit vectors in the directions of r and 0, respectively. (See 
fig. 1.) 

The differential equation that defines the streamline geometry may now be formed 
from equation (9) as 
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9 <f>\ _ • ^0 _ -(1 + r , 2)sin 0 sin 


\9x g g _ j cos o!+ r' cos 0 sin a 

Although the left side of this equation is written as a partial derivative, (3 is constant 
along a streamline and therefore it may be replaced by an ordinary derivative for the 
integration of a given streamline. Thus, for any chosen streamline emanating from the 
stagnation point, equation (11) may be integrated to determine the geometry of this stream- 
line (that is, 0 = 0(x)) as it wraps around the body. Note that the modified Newtonian 
pressure distribution locates the stagnation point at r 0 f = cot a, 0 = it radians. 

Scale Factor 

The next task is to determine a relation for the scale factor h along a streamline. 
Let S |3 be a unit vector in the direction of /3; thus 


e /3 - e s x e n 


(i + e r r')sin a sin 0 + e^(r' sin a cos <p + cos a) 


Now since the element of arc length in the /3 direction on the surface is h d/3, then 


h/ 8 j3\ _ ~ - _ r* sin a cos 0 + cos a 

= e f> ' * ‘ A 


Therefore, along the surface 


r(r' sin a cos 0 + cos a) / 90' 


where 0 = 0(x,/ 3) defines the circumferential location of the inviscid surface stream- 
lines, and again /3 is constant along a given streamline. 

Since (p = cp{x,j3) is analytic, 

9^0 920 

9x 9/3 ~ 9/3 9x 

” d since $) x ■ (4) j$) x ’ 


9 2 0 _ p 9 / 9 0\1 ( 9 0 


9/3 9 x ^0\9xyy \ 9 ^/ x 
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Then equation (15) becomes 


Alnl^l 


8 ( d 4 >\ 

3x \BB I 
X K/ X 

£ 

0 4>\ 9 x/ |3 


x 


( 17 ) 


Taking the partial derivative of equation (11) with respect to $ and holding x constant 
yields 

-(1 + r'2)(r' tan a + cos <£)tan a 
r(l + r T cos <p tan a ) 2 



With this substitution, equation (17) finally becomes 



-(1 + r'2)(r T tan a + cos <p ) tan a 
r(l + r* cos <p tan a) 2 


(18) 


For the integration of equation (18) along a given streamline, the left side of this equation 
may be written as an ordinary derivative: 


d_ 

dx 



-(1 + r'2)(r’ tan a + cos <£)tan a 
r(l + r' cos <p tan a) 2 


(19) 


The solution of this equation allows h to be computed from equation (14). 

Equations (2) and (4) may also be used to compute the pressure ratio and velocity 
ratio along any inviscid surface streamline for the body at an angle of attack. For equa- 
tion (2), cos i// may be evaluated from equations (7) and (8) as 


cos \ p = 


V„ 


e n 


r' cos of - cos <p sin a 
(1 + r’ 2 ) 1 / 2 


(20) 


However, the surface region on the leeward side of the body at an angle of attack where 
cos \j/ = 0 will be in the "aerodynamic shadow," and the modified Newtonian theory does 
not predict the pressures in this region. In the present method the surface pressures in 
the shadowed region are assumed to be free-stream static pressure, but the streamline 
geometry and the scale factor h are still determined from equations (11), (14), and (19) 
in this region. 


Application to Spherically Blunted Cone at Angle of Attack 

The determination of the streamline geometry and the scale factor h is very sim- 
ple for spherically capped bodies at an angle of attack when the surface sonic line (invis- 
cid) lies entirely on the spherical surface. In this case the streamlines on a portion of 
the spherical cap follow spherical meridians about an axis through the center of the sphere 
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and parallel to V^. From geometrical considerations, the coordinate /3 on the spherical 
cap is the angular position of a streamline measured in the plane tangent to the spherical 
cap at the stagnation point, with /3 = 0 for the top meridian and 0 = ir radians for the 
bottom meridian. (See fig. 2.) Also, while on the spherical cap, the scale factor h is 
simply the distance from the point in question on the sphere to the axis through the cen- 
ter of the sphere and parallel to V^. Therefore, for the spherical cap, 

1 1/2 


— = sin xl/ = 

R 


r'2 + 1 - (r* cos a - cos <b sin a)2 


<2 + 1 


( 21 ) 


and from the geometry of the streamlines, it may be verified that at the sphere-cone 
interface /3 and <p^ are related by the expression 


Hence it follows that 

cos /3 = 


hi sin /3 = rj sin 

cos <£i + k tan a 


(1 + k2)sin2<£i tan^o!+ (l + k cos (pi tan 


( 22 ) 


where i refers to values at the sphere-cone interface and k is the slope of the cone 
surface: 

k = tan 9 = r’ (23) 

For the conical afterbody, equation (11) may be integrated to yield the streamline 
geometry as , 


r. 



\ (l+k2)tan a 


fsin (p-\ 
^sin <p j 


k2 

l+k2 


(24) 


To determine the scale factor h for the afterbody, it is noted that 


80 

9/3 


/ d(j) \ /dtp A 

” (a^Tj \bb~) 

X \ r V x \ / x 


(25) 


and taking the partial derivative of equation (24) with respect to 4> ^ (holding x constant) 
yields 


' d(p\ A + k tan a cos S in 
90i/ x 1 1 + k tan a cos cp Jsin <p^ 


(26) 


Now equation (22) is differentiated with respect to /3 (holding x constant) to get 
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sin2^ i + (cos cos a + k sin a^‘ 
cos a + k sin a cos 


( 27 ) 


Finally, for the conical afterbody, equations (24) to (27) are substituted into equation (14) to 
obtain 

k 

/ 4>\( l+k 2 )tano! o 

h / tan 2| /sin c/A 1+k sin20 i + (cos cos a + k sin aj 

R ytan \ n V [sin2<£ + (cos 4> cos a + k sin + k^) 1 ^ 

(28) 

The reduction of equation (28) for the limiting condition of a = 0 and for the case of 
<p = 0, xr with a finite is given in the appendix. 

Equations (24) and (28) may now be used in equations (1) to (4), with r replaced by 
h and s measured along an inviscid surface streamline, to determine the heating-rate 
ratio along any streamline. The integral in the denominator of equation (1) can be evalu- 
ated by any simple technique, such as Simpson's rale. 


RESULTS 


The method presented here was used to compute the nondimensional heating rate 
over spherically blunted cones with 15° and 30° semiapex angles. Angles of attack of 0°, 
10°, and 20°, a Mach number of 10.6, and a ratio of specific heats of 1.4 were considered. 
Calculated variations have been compared with the experimental data given in reference 6. 

Figures 3 and 4 give the longitudinal variation of the nondimensional heating rate 
along the bottom, top, and side (<p = radians) at angles of attack of 0°, 10°, and 20°. It 
is clear from these figures that theoretical and experimental data agree well along the 
windward and side meridians, but along the top meridian (leeward side) the agreement is 
poor. Also, near the spherical cap theoretical values generally exceed experimental data. 
This would seem to indicate that the spherical-cap heating-rate distributions are not pre- 
dicted accurately near the sphere- cone interface. 

Variations around the bodies at specific longitudinal locations are given in figures 5 
and 6. Previous observations made concerning the longitudinal variations are again valid; 
that is, the agreement between experimental and theoretical data is worst on the leeward 
side and near the nose cap. 

It should be noted that when the angle of attack is large enough to cause any portion 
of the sonic line to move off the spherical cap, the stagnation point will move a short dis- 
tance from the position r 0 ' = cot a, <p = n radians (see ref. 7). In addition, the axial 
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symmetry on a portion of the spherical cap, about an axis through the center of the sphere 
and parallel to V , is destroyed. The present method does not account for this phenome- 
non, and the spherically blunted 30° cone at a = 20° is a marginal case. 


CONCLUDING REMARKS 

A simplified method for computing the laminar heating rates over bodies at an angle 
of attack is developed and the results are shown to agree well with experimental results 
for the case of a spherically blunted cone. Although the applications presented here are 
blunted cones, the basic method is applicable to any asymmetric or axisymmetric body at 
an angle of attack. The relatively small amount of numerical computations required 
makes the method particularly attractive for engineering applications. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., June 7, 1968, 

124-07-01-18-23. 
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APPENDIX 


LIMITING FORM OF SCALE FACTOR h FOR a = 0 
AND FOR (j) = 0,ir WITH a > 0 

Equation (24) may be used to rewrite equation (28) in the following form: 
h r sin cj) sin2$ ^ + ( cos <f>^ cos a + k sin a) 1 


R i sin cp j + ( cos 0 cos a + k sin a)2j^^(l + k2)*/^ 

Then if equation (24) is rearranged to give 

il± k2 > tan a tan 0i /- ^\ ktan « 
k = tan Tf sm 0i\ 

’ tan — \ sin v 
2 


it is clear that in the limit as a — 0, 


1 = 


. 0 i 

tan — 

4 


or 


0 i = 0 

With this identity substituted into equation (Al), there results 

f h' 


(Al) 


(A2) 


(A3) 


( R o? —0 r i(l + k2)l/2 

But r^ = R(1 + k2) - ^/^ and the correct limit for h as a — 0 is obviously obtained: 

h G-=0 = r 

The next problem is to find the limit of h as <p — 0 and as p — u radians (for 
a. > 0). From equation (Al), it is found for the first limit that 


l im (U\ = JL (1 + k tan ojcos a /sin </>' 


Aim — ~ — urn | 

r ‘ (l + k2) 1/2 f-oA sln ^v 

0i~O 0 i-0 


(A4) 
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APPENDIX 


In order to evaluate lim [ sin . ft V the trigonometric identity tan % = — S * n — — is sub- 

<p-0 \sin (pi J 2 i + cos <p 

cPi-0 


stituted into equation (24) to get 


k 


k+k2tan ct 


r_ = / 1 ± cos 0\ (1+k2)tan “/sin_^\(l+k 2 )tan a 
r i yl + cos \sin (p ) 


(A5) 


Then from equation (A 5), 


lim 


sin (p 


cp - 0 \ sin <t > i 

4>i~0 



(l+k2)tan a 
k+k2tan a 


Now substitute this relationship into equation (A4) to get 

k-tan a 

lim /h \ - / r \k+k2tan a (l + k tan a) cos a 

\ r i / (1 + k 2) X / 2 

(Pi-0 

Next consider lim (— •). From equation (Al), 

(p- 7r' R/ 

(Pi-Tl 


lira (!) = i (1 - fe tan deoe a Um /sta£\ 
(P-tjKR/ r i (l + k2)-*-/ 2 it \ s * n $t) 

(pi~V (pi-TT 


(A6) 


(A7) 


Note that since the stagnation point must lie on the spherical cap and since r 0 ' = cot a 
at the stagnation point, k tan a < 1 everywhere on the conical afterbody. 


Now the trigonometric identity 

to get 



- — ^° S — is substituted into equation (24) 
sin <p 



( 1 - cos <ftj\ 

\1 - COS (p ) 


k 


(l+k2)tan a 


fein 0A 
k sin (p ) 


k2tan a-k 
(l+k2)tan a 
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Then in the limit this equation goes to 


(l+k2)tan a 


lim /sin <A J r\ k-k 2 tan a 
0— it \ sin 4>iJ \ r i / 


(A8) 


Finally, equation (A8) is substituted into equation (A7) to obtain 


k+tan a 


lim (r) = (n) 
p-Tr\ R / \ r v 


k-k^tan a (i _ k tan o') cos a 
1 (1 + k2)V2 


Therefore, equations (A6) and (A9) give the limiting values of at <p = 0 and <fi = ir 

*v 

radians, respectively. It should be noted that equation (A9) gives the same relation as 
equation (A6) when a is replaced with -a. 
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Figure 2.- Sketch illustrating streamline coordinate p. 
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Figure 3.- Longitudinal variation of nondimensional heating- rate ratio q w ,o on a spherically blunted cone with 15° semiapex angle. IVU 10.6. 
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